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1 Introduction 



Consider a sequence of random variables {Fn,n > 1} defined on a complete 
probability space {0,,J^,P). Suppose that the cj-field T is generated by an 
isonormal Gaussian process X = {X{h), h G f)} on a real separable infinite- 
dimensional Hilbert space S^. This just means that X is a centered Gaussian 
family of random variables indexed by the elements of Sj, and such that, for 
every h,g G S), 

E[X{h)X{g)]={h,g)s,. (1.1) 

Suppose that the sequence {Fn,n > 1} is normalized, that is, E{Fn) = 
and lim„_+oo -£'(-^n) = 1- A natural problem is to find suitable conditions 
ensuring that F„ converges in law towards a given distribution. When the 
random variables Fn belong to the qth Wiener chaos of X (for a fixed q >2), 
then it turns out that the following conditions are equivalent: 

(i) Fn converges in law to A^(0, 1); 
'41 - 3; 



(ii) lim„_oo-E^[^„ 



(iii) lim„_^oo \\DF, 



q in L^{n). 



Here, D stands for the derivative operator in the sense of Malliavin calculus 
(see Section [2] below for more details). More precisely, the following bound 
is in order, where denotes a standard Gaussian random variable: 



sup\P{Fn < z) - P{N < z)\ ^ 



E 



1 



l--\\DFn\\l 



(1.2) 



^ J^VI^ra-3|- (1-3) 



The equivalence between conditions (i) and (ii) was proved in Nualart 
and Peccati [22] by means of the Dambis, Dubins and Schwarz theorem. 
It implies that the convergence in distribution of a sequence of multiple 
stochastic integrals towards a Gaussian random variable is completely de- 
termined by the asymptotic behavior of their second and fourth moments, 
which represents a drastic simplification of the classical "method of mo- 
ments and diagrams" (see, for instance, the survey by Peccati and Taqqu 
[26], as well as the references therein). The equivalence with condition (iii) 
was proved later by Nualart and Ortiz-Latorre [21] using tools of Malliavin 
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calculus. Finally, the Berry-Esseen's type bound (|1.2p is taken from Nour- 
din and Peccati |16j . while (jl.Sp was shown in Nourdin, Peccati and Reinert 

Peccati and Tudor [27] also obtained a multidimensional version of the 
equivalence between (i) and (ii). In particular, they proved that, given a 
sequence {Fn,n > 1} of d-dimensional random vectors such that belongs 
to the Qith Wiener chaos for i = 1,. . . ,d, where 1 ^ qi ^ . . . ^ q^, then 
if the covariance matrix of -F„ converges to the d x d identity matrix I^, 
the convergence in distribution to each component towards the law A^(0, 1) 
implies the convergence in distribution of the whole sequence F„ towards 
the standard centered Gaussian law A^(0, I^). 

Recent examples of application of these results are, among others, the 
study of p- variations of fractional stochastic integrals (Corcuera et al. [3]), 
quadratic functionals of bivariate Gaussian processes (Deheuvels et al. [5] ) , 
self-intersection local times of fractional Brownian motion (Hu and Nu- 
alart [7]), approximation schemes for scalar fractional differential equa- 
tions (Neuenkirch and Nourdin [12] ) , high-frequency CLTs for random fields 
on homogeneous spaces (Marinucci and Peccati [10^ [IT] and Peccati [23]), 
needlets analysis on the sphere (Baldi et al. [1] ) , estimation of self-similarity 
orders (Tudor and Viens [H]), weighted power variations of iterated Brow- 
nian motion (Nourdin and Peccati [15]) or bipower variations of Gaussian 
processes with stationary increments (Barndorff-Nielsen et al. [2]). 

Since the works by Nualart and Peccati ^22j and Peccati and Tudor [2^ , 
great efforts have been made to find similar statements in the case where 
the limit is not necessarily Gaussian. In the references }24j and [25], Peccati 
and Taqqu propose sufficient conditions ensuring that a given sequence of 
multiple Wiener- Ito integrals converges stably towards mixtures of Gaussian 
random variables. In another direction, Nourdin and Peccati [Tl] proved 
an extension of the above equivalence (i) - (iii) for a sequence of random 
variables {Fn,n > 1} in a fixed qth Wiener chaos, q > 2, where the limit 
law is 2Gyi2 — 1^, being the Gamma distribution with parameter z^/2. 

The purpose of the present paper is to study the convergence in distri- 
bution of a sequence of random variables of the form Fn = S'^{un), where 
Un are random variables with values in Sj^'^ (the qth tensor product of S^) 
and 5'^ denotes the multiple Skorohod integral (that is, 6'^{u) = d{5{u)), 
(5^(n) = 6{5{S{u))), and so on), towards a mixture of Gaussian random vari- 
ables. Our main abstract result. Theorem 13. 11 roughly says that under some 
technical conditions, if {un, D'^Fn)^0q converges in L^(n) to a nonnegative 
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random variable S'^, then the sequence Fn converges stably to a random vari- 
able F with conditional characteristic function E (^e^^^^ X) = E{e~~^ ). 
Notice that if Un is deterministic, then Fn belongs to the gth Wiener chaos, 
and we have a sequence of the type considered above. In particular, if 5^ is 
also deterministic, we recover the fact that condition (iii) above implies the 
convergence in distribution to the law A^(0, 1). 

We develop some particular applications of Theorem 13. II in the following 
directions. First, we consider a sequence of random variables in a fixed 
Wiener chaos and we derive new criteria for the convergence to a mixture 
of Gaussian laws. Second, we show the convergence in law of the sequence 
where q>2 and Un is a g-parameter process of the form 

n-l 

Un = f{Bk/n)'^ik/n,(k+l)/n]i ; 

A:=0 

towards the random variable aH,q Jq f{Bs)dWs, where B is a. fractional 
Brownian motion with Hurst parameter H £ ^ is a standard 

Brownian motion independent of B, and crH,q denotes some positive con- 
stant. This convergence allows us to establish a new asymptotic result for 
the behavior of the weighted gth Hermite variation of the fractional Brow- 
nian motion with Hurst parameter H E ^^j^^! which complements and 
provides a new perspective to the results proved by Nourdin [13], Nour- 
din, Nualart and Tudor [18], and Nourdin and Reveillac [19j . The reader is 
referred to Section [5] for a detailed description of these results. 

The paper is organized as follows. In Section [21 we present some prelimi- 
nary results about Malliavin calculus. Section [3] contains the statement and 
the proof of the main abstract result. In Section [U we apply it to sequences 
of multiple stochastic integrals, while Section [5] focuses on the applications 
to the weighted Hermite variations of the fractional Brownian motion. 

2 Preliminaries 

Let ^ be a real separable infinite-dimensional Hilbert space. For any integer 
g > 1, let ij®* be the gth tensor product of ^. Also, we denote by i^®'' the 
gth symmetric tensor product. 

Suppose that X = {X{h), h £ Sj} is an isonormal Gaussian process on 
Sj, defined on some probability space {0,,T,P). Recall that this means that 
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the covariance of X is given in terms of the scalar product of Sj by (jl.ip . 
Assume from now on that is generated by X. 

For every integer q > 1, let Tig be the qth Wiener chaos of X, that 
is, the closed linear subspace of L'^{Q) generated by the random variables 
{Hq{X{h)),h G -^3)11^11^ = 1}, where Hq is the gth Hermite polynomial 
defined by 

We denote by TLq the space of constant random variables. For any g > 1, 
the mapping Iq{h'^'^) = q\Hq(X{h)) provides a linear isometry between ^®'' 
(equipped with the modified norm \/^ Ihllj^®?) ^^'^ '^q (equipped with the 
L^(f2) norm). For g = 0, by convention TLq = R, and Iq is the identity map. 

It is well-known (Wiener chaos expansion) that L^(r2) can be decom- 
posed into the infinite orthogonal sum of the spaces TLq. That is, any square 
integrable random variable F G Lp'iQ.) admits the following chaotic expan- 
sion: 

oo 

i^ = E^'?(/^)' (2-1) 

where /o = E[F], and the fq G ^®^, > 1, are uniquely determined by 
F. For every g > 0, we denote by Jq the orthogonal projection operator 
on the qih. Wiener chaos. In particular, if F G L^(r2) is as in (12. ip . then 
JqF = Iqifq) foT cvcry q > 0. 

Let {cfc, fc > 1} be a complete orthonormal system in S). Given / G ^®^, 
g G ij®*^ and r G {0, . . . ,pAq}, the rth contraction of / and g is the element 
of i^^Cp+'y-Sr) defined by 

oo 

f®rg= E (/, (g) . . . (g) ei,,)^®r (g) (5,eii (g) . . . (g) ei,,)i5®r. (2.2) 



Notice that f ®r g is not necessarily symmetric. We denote its symmetriza- 
tion by f®rg £ ^©(p+'J-^r)^ Moreover, f ®o g = f ® g equals the tensor 
product of / and g while, iov p = q, f (dq g = (/, g)sj<si- 

In the particular case $) = L'^{A,A,fi), where {A, A) is a measurable 
space and ^ is a u-finite and non-atomic measure, one has that ^®'? = 
Lg{A'^ , A'^'^ , fi'^'^) is the space of symmetric and square integrable functions 
on A'l. Moreover, for every / G i^®"^, Iq{f) coincides with the multiple 
Wiener-Ito integral of order q of f with respect to X (introduced by Ito in 
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[5]) and (|2.2p can be written as 



if '^r g){tl, . . . ,tp+q-2r) = / f {ti, . . . ,tp-r , Si, . . . , Sr) 

X g{tp_r+l, . . . ,tp+g_2r,Sl, ■ ■ . , Sr)dn{si) . ..dfJ,{Sr). 

Let us now introduce some basic elements of the Malliavin calculus with 
respect to the isonormal Gaussian process X. We refer the reader to Nualart 
|20| for a more detailed presentation of these notions. Let S be the set of 
all smooth and cylindrical random variables of the form 

F = g{X{4>i),...,X{4>n)), (2.3) 

where n > 1, g : M" ^ M is a infinitely differentiable function with compact 
support, and (pi £ f). The Malliavin derivative of F with respect to X is the 
element of L^(J7,i5) defined as 

1=1 * 

By iteration, one can define the qth derivative D''F for every q > 2, which 
is an element of 

For q > 1 and p > 1, P"^'^ denotes the closure of S with respect to the 
norm || • Wm-p, defined by the relation 



IP 

\ni,p 



e[\f\p] + Y,e(\\d'f 



i=l 



The Malliavin derivative D verifies the following chain rule. If ip : 

is continuously differentiable with bounded partial derivatives and if F 

{Fi, . . . , Fn) is a vector of elements of B-*^'^, then ^{F) G ID^'^ and 

dip 



1=1 



We denote by 6 the adjoint of the operator D, also called the divergence 
operator. The operator 6 is also called the Skorohod integral because in the 
case of the Brownian motion it coincides with the anticipating stochastic 
integral introduced by Skorohod in [30J. A random element u S L^(J7,^) 
belongs to the domain of 5, noted T)om5, if and only if it verifies 



\E{{DF,u)^) \ < Cu VE{F^) 
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for any F G D^'^, where is a constant depending only on u. If u G Domd, 
then the random variable 5{u) is defined by the duality relationship (called 
'integration by parts formula'): 

E{F5{u))=E{{DF,u)^), (2.4) 

which holds for every F G B^'^. The formula ()2.4p extends to the multiple 
Skorohod integral 5^, and we have 

E {F5'i{u)) = E {{DIE, (2.5) 

for any element u in the domain of 5^ and any random variable F G B'^'^. 
Moreover, 6'^{h) = Iq{h) for any h G ^®^. 

The following property will be extensively used in the paper. 

Lemma 2.1 Let q > 1 be an integer. Suppose that F G D'^'^, and let u 
be a symmetric element in Dom6'' . Assume that, for any < r + j < 
q, {D-F,6^{u))^^,. G L2(17,^®'?— i). Then, for any r = 0, . . . ,q - 1, 
{D^'F,u)^^r belongs to the domain of 5''^^ and we have 

FS^iu) = ^2 ('^) S'"'' {{D^F, u)^^r) . (2.6) 

(We use the convention that 5^{v) = t; G M, and D^F = F, F £ L^(0).j 

Proof. We prove this lemma by induction on q. For g = 1 it reads 
F5{u) = 5{Fu) + {DF,u)sj, and this formula is well-known, see e.g. [201 
Proposition 1.3.3]. Suppose the result is true for q. Then, if u belongs to 
the domain of (5'^^^, by the induction hypothesis applied to 5{u), 

F5^+i(n) = FS^idiu)) = j2 f^)^'"'' {{D^F,6iu))^^,.) . (2.7) 

On the other hand 

{D'E, 6iu))^^r = S {{D'E, n)^«.) + {D'+^E, u) . (2.8) 

Finally, substituting (12. Sp into (12. 7p yields the desired result. ■ 

For any Hilbert space V, we denote by I])^'P{V) the corresponding Sobolev 
space of F- valued random variables (see ^(H page 31]). The operator 6*^ 
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is continuous from I])^'P(^S^^'^) to D'^ '''^j for any p > 1 and any integers 
k > q > 1, that is, we have 

P''(^)IId'=-'?.P — '^k,p\\u\\jj,k,p(^^0q-^ (2-9) 

for ah u G I}^''p[S)^'^), and some constant c^^p > 0. These estimates are con- 
sequences of Meyer inequaUties (see j20l Proposition 1.5.7]). In particular, 
these estimates imply that B'^'^(i3®'') C DomJ"? for any integer q > 1. 

We will also use the following commutation relationship between the 
Malliavin derivative and the Skorohod integral (see \20\ Proposition 1.3.2]) 

D6{u) = u + 6{Du), (2.10) 

for any u G B^'^(^). By induction we can show the following formula for 
any symmetric element u in W^^''^{Sj^^) 

D^5\u) = ^ (^^iW-\D^-'u). (2.11) 

We will make use of the following formula for the variance of a multiple Sko- 
rohod integral. Let u,v G B^'''^(i3'^'') C Domd'^ be two symmetric functions. 
Then 

1 / \ 2 



E(') ^lE{{n,6''-\D'>'^v))^, 
j2 (fj '^lE [{D^-'u, D''~^'v)^^,,,_^,) . (2.12) 

i—n V / 



i=0 

The operator L is defined on the Wiener chaos expansion as 

oo 

L = ^-qJ^, 

and is called the infinitesimal generator of the Ornstein-Uhlenbeck semi- 
group. The domain of this operator in L'^{Q,) is the set 



DomL = {F e L^{VL) -.^q^ Pq^Wl^^^^) < oo} 

q=l 



^2,2 
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There is an important relation between the operators D, 5 and L (see |20l 
Proposition 1.4.3]). A random variable F belongs to the domain of L if and 
only if F e Dom {6D) (i.e. F G ID^'^ and DF G Domd), and in this case 

5DF = -LF. (2.13) 

Note also that a random variable F as in (|2.1|) is in B^'^ if and only if 

oo 
q=l 

and, in this case, E {\\DF\\l) = Eg>iQQ^-\\fq\\l^,- If ^ = L^{A,A, fi) (with 
fj, non-atomic), then the derivative of a random variable F as in (j2.ip can 
be identified with the element of L'^{A x Cl) given by 

oo 

DaF = Y,qIg^i{f,{;a)), a e A. (2.14) 

9=1 

Finally, we need the definition of stable convergence (see, for instance, 
the original paper [29], or the book [9] for an exhaustive discussion of stable 
convergence) . 

Definition 2.2 Let F^ be a sequence of random variables defined on the 
probability space {^},J-,P), and suppose that F is a random variable defined 
on an enlarged probability space {0,,Q,P), with Q Q. We say that F^ 
converges ^-stably to F (or only stably when the context is clear) if, for any 
continuous and bounded function / : M ^ M and any bounded T -measurable 
random variable Z , we have E [f{Fn)Z\ E [f{F)Z\ as n tends to infinity. 

3 Convergence in law of multiple Skorohod inte- 
grals 

As in the previous section, X = {X(h),h G f)} is an isonormal Gaussian 
process associated with a real separable infinite-dimensional Hilbert space 
S^. The next theorem is the main abstract result of the present paper. 

Theorem 3.1 Fix an integer q > 1, and suppose that Fn is a sequence of 
random variables of the form Fn = 5'^{un), for some symmetric functions Un 
in ©^'^'^'^(i^'^''). Suppose moreover that the sequence Fn is bounded in L^{il), 
and that: 
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(i) {DFn f^' ... (Di-^Fn)^''''-' ® /i>^«, converges in L\Q) to 
zero, for all integers r,ki, . . . , kq-i > such that 

ki+2k2 + ... + {q- l)kq-i +r = q, 

and all h G 

(ii) {un, W^Fn) ^cgiq converges in L^(Q) to a nonnegative random variable 

Then, Fn converges stably to a random variable with conditional Gaus- 
sian law N{0, S'^) given X. 

Remark 3.2 When q = 1, condition (i) of the theorem is that {un,h)sj 
converges to zero in L^(r2), for each h £ S). When q = 2, condition (i) means 
that {un-,h® g) fjfs2, {un,DFn'Sih)^02 and {un, DFn DFn) converge to 
zero in L^(Q,), for each h,g S Sj. And so on. 

Proof of Theorem 13.11 Taking into account Definition 12.21 it suffices to 
show that for any hi, ... , hm G S^, the sequence 

in = {Fn,X{hi),...,X{hra)) 

converges in distribution to a vector [Fao,X{hi), . . . ,X{hm))-, where Fqo 
satisfies, for any A £ M, 

E{e^^^-'\X{hi),...,X{hm))=e-'^'\ (3.1) 

Since the sequence F„ is bounded in the sequence is tight. As- 

sume that {Foo,X{hi), . . . ,X[hm)) denotes the hmit in law of a certain 
subsequence of denoted again by ^n- 

Let Y = (t){X{hi),...,X{hm)), with 4) G C^{W^) (0 is infinitely dif- 
ferentiable, bounded, with bounded partial derivatives of all orders), and 
consider 0„(A) = E {e'^^^Y) for A e M. The convergence in law of , 
together with the fact that Fn is bounded in L^(r2), imply that 

lim (/)'„(A) = lim iE ( Fne'^^'^Y^ = iE{F^e'^^^Y). (3.2) 

n— »oo n— >oo V / 
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On the other hand, by ()2.5p and the Leibnitz rule for D'^, we obtain 
= iE((un^D^{e'^^-Y] 



a=0 

xE (e'^^- (^un, {DFnf'''^ . . . ®{D''FnT^-'®D'''''Y^ ^ J 



xE fe'^-^" / n„, iDFnf^' ... (D'^Fn)®^- D^-^^y 

V \ / 

where the second sum in the two last equalities runs over all sequences of 
integers (fci, . . . , kg) such that ki + 2k2 + . . . + aka = a, due to the Faa di 
Bruno's formula. By condition (i), this yields that 

ct>'^{\) = -XE (e*^^" {un, D'iFn)^^, y) + i?„, 

with Rn converging to zero as n — > cxd. Using condition (ii) and (|3.2|) . we 
obtain that 

iE{Fooe'^^^Y) = -XE [e'^^°° S'^Y^ . 

Since S"^ is defined through condition (ii), it is in particular measurable with 
respect to X. Thus, the following linear differential equation verified by the 
conditional characteristic function of F^o holds: 

^E{e'^^^\X{hi), X{hm)) = -X E{e'^^-'\X{hi), X{h^)). 



By solving it, we obtain ()3.ip . which yields the desired conclusion. ■ 

The next corollary provides stronger but easier conditions for the stable 
convergence. 

Corollary 3.3 For a fixed ^ > 1, suppose that Fn is a sequence of random 
variables of the form Fn = 5'^{un), for some symmetric functions Un in 
]I])2g,2g^^®g-j^ guppggg moreovcr that the sequence Fn is hounded in W^'p for 
all p>2, and that: 
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(i') {un, h)fji)t)q converges to zero in L^{VL) for all h G S^®'^ ; and u„ Cg)/ D^F^ 
converges to zero in ^®('^~')) for all I = 1, . . . ,q — 1; 

(ii) {un, D'^Fn) converges in L^{Q) to a nonnegative random variable 

Then, Fn converges stably to a random variable with conditional Gaus- 
sian law N{0, S'^) given X. 

Proof. It suffices to show that condition (i') imphes condition (i) in 
Theorem 13. II When /c^ 7^ for 1 < a < — 1, we have, for aU h E (with 
r = q - ki - 2k2 - . ■ ■ - aka), 



Un ®a D^Fn, 



I f|8>9 



® (D'^-'^Fr,)®^--^ ® {D'^Fnf^^'^-^^ ® h 



< WUr, 



D^F 



n\\f}0(.q-a) 



{DFnf"' 



{D^-^Fn)®^^'^ ® {D''Fnf^^^~^^ ® h 



i(q-a) 



f^^(q-a) 



The second factor is bounded in L^(r2), and the first factor converges to zero 
in L^(il), for all a = 1, . . . , g — 1. In the case a = we have that (u^, 
converges to zero in L^(Q,), for all h £ S)®'^, by condition (i'). This completes 
the proof. ■ 



4 Multiple stochastic integrals 

Suppose that is a Hilbert space L'^{A,A,fi), where {A, A) is a measurable 
space and /i is a c-finite and non-atomic measure. 

Fix an integer m > 2, and consider a sequence of multiple stochastic 
integrals {F^ = Im{gn), n > 1} with g^ G i^®"^. We would like to apply 
Theorem 13.11 with g = 1 to the sequence Fn- To do this, we represent each 
Fn as 

Fn = 5{Un), with Un = Im-lidn), 

for gn € Sj®^ some function which is symmetric in the first m — 1 variables. 
Notice that, from (j2.14p . we have DFn = mlm~iign)- Hence, since 
= ~m^^n = m^^^^-n) (EZS]), gn IS always a possible choice for gn- 

(In this case, gn is symmetric in all the variables.) However, as observed. 
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for instance, in Example 14.21 below, the choice gn = gn does not allow to 
conclude in general. 

Proposition 4.1 For a fixed integer m > 2, let Fn he a sequence of random 
variables of the form F„ = Imign), with gn S Suppose moreover that 

Fn is bounded in L'^{Q) and that Fn = S{un), where Un = Im-iign); for 
dn £ i^®™ some function which is symmetric in the first m — 1 variables. 
Finally, assume that: 

(o-) {dn ®m-i dn, h'^'^) ^<g,2 convcrgcs to zero for all h £ S^; 

(b) {un, DFn) sj converges in L^(Q) to a non negative random variable S*^. 

Then, Fn converges stably to a random variable with conditional Gaus- 
sian law N{0, 5^) given X. 

Proof. It suffices to apply Theorem 13.11 to Un = Im-iidn) and q = 1- 
Indeed, we have 

E{{Un,h)l) = =S(/„_i(5„®i/l)2) 

= (m - 1)! (g)i 

= (m - 1)! (5„ 5„,/i®^)^552 ^ 0, 

which implies condition (i) in Theorem 13. H see also Remark 13. 2[ Condition 
(ii) in Theorem 13.11 follows from (b). ■ 

Example 4.2 (see also [28] Proposition 2.1] or [24^ Proposition 18] for two 
different proofs using other techniques). Suppose that {Wt, t S [0,1]} is 
a standard Brownian motion. (This corresponds to A = [0, 1] and ^ the 
Lebesgue measure.) Assume that m = 2 and take gn{s,t) = ^^/n{s V t)^. 
Then ^ 

Fn = h{9n) = V^ [ t^WtdWu 
Jo 

and ^ 

DsFn = V^s'^Ws + t^WtdWt. 

We can take Un{t) = ^/nt'^Wt, that is, gnis,t) = ^/nt^l[Q^t]{s). In this case, 

{dn (S)ign){s,t) = ns''r(s At), 
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which converges to zero weakly in L^(r2), and 



{Un,DFn)sj 



s'^W.dW, dt 



which converges in Lp'[^) to ^Wf. Therefore, conditions (a) and (b) of 
Proposition 14.11 are satisfied with = ^W^, and F„ converges in distri- 
bution to ^M^i X iV, with iV ~ Af(0, 1). One easily see on this particular 
example that the choice gn = gn does not allows us to conclude in gen- 
eral (except when 5^ is deterministic); indeed, one can check here that 
{un, DFn)f} = ^\\DFn\\j^ docs uot Converge in L^(0). 

If we take gn = gn and = 1, then condition (b) coincides with con- 
dition (iii) in the introduction. In this case, Nualart and Peccati criterion 
combined with Lemma 6 in [21] tells us that, if the sequence of variances 
converges to one, then condition (a) is automatically satisfied. 

On the other hand, we can also apply Theorem 13.11 with m„ = gn- In this 
way, applying Corollary 13.31 we obtain that the following conditions imply 
that Fn converges to a normal random variable A^(0, 1) independent of X: 

(a) gn converges weakly to zero; 

(/?) Wdn^i gn\\fj<s2{q-t) converges to zero for alH = 1, . . . , g — 1; 

(7) ^'Ibnll^®, converges to 1. 

Indeed, notice first that if gn is bounded in S)®'^, then Fn is bounded 
in all the Sobolev spaces D'^'^, p > 2. Then, condition (ii) in Corollary 13.31 
follows from (7) and the equality {Iq{gn)) = q\gn- On the other hand, 
condition (i') in Corollary 13.31 follows from (ii) and 



E 



D'Fn 



2 

f, '8(9-0 



7I2 



{q-iy. 
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|2 



(q-iy? 

nl2 



E \\gn<^l Iq-l{gn)\\^s>{q-l) 
E \\Iq-l{gn®l gn)\\\is>(q-l) 



(^q — \\9n^ign\\^(S>2(q-l) 



7l2 



< 



_ \\gn ®i gn\\ ^<S>2(q-l) ■ 



In this way we recover the fact that condition (iii) in the introduction implies 
the normal convergence. 
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5 Weighted Hermite variations of the fractional 
Brownian motion 



5.1 Description of the results 

The fractional Brownian motion (fBm) with Hurst parameter H G (0, 1) is 
a centered Gaussian process B = {Bt, t >0} with the covariance function 

E{BsBt) = Rh(.s, t) = ^ (t^^^ + -\t- . (5.1) 

From dSH), it follows that E\Bt - Bs\^ = {t - sf^ for all ^ s < t 
and that, for each a > 0, the process {a~^ Bat, i > 0} is also a fBm with 
Hurst parameter H (self-similarity property). As a consequence, the se- 
quence {Bj — Bj-i, j = 1,2, . . .} is stationary, Gaussian and ergodic, with 
correlation given by 

p^in) = i [In + Ip^ - 2|np^ + In - ll^^] , (5.2) 

which behaves as H{2H — l)|np^~^ as n tends to infnity. 

Set ABk/n = S(fc+i)/n - ^fc/n, where /c = 0, 1, . . . , n, and n > 1. The 
ergodic theorem combined with the self-similarity property implies that the 
sequence n'^^~^ J^kZo (^-Sfc/n) converges, almost surely and in L^(Q), to 
E{Bf) = 1. Moreover, it is well-known (see, e.g., [3]) that, provided H G 
(0, |), a central limit theorem holds: the sequence 

^ n—l ^ n—1 

-j= (n^^ {AB,/^Y - l) = 4 E ^2 (n^A5,/„) (5.3) 
^ k=0 ^ fc=0 

converges in law to N{0,a'jj) as n ^ oo, for some constant an > 0. (Notice 
also that, by normalizing with y'nTogn instead of y/n, the central limit 
theorem continues to hold in the critical case H = |.) When F > f , the 
situation is very different. Indeed, we have in contrast that 



n—l n—l 



{AB,/^f - l) = ni-2^ Yl ^2 (n^AB,/J 



k=0 k=0 



converges in L'^{Q). More generally, consider an integer q > 2. If H < 1 — 
then the sequence 

— ^ (n^Ai?,/„) (5.4) 
^ k=0 
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converges m law to A^(0,cj2^) (for some constant CTq^H > 0), whereas, if 
H > 1 — then the sequence 



ra-l 



n 

k=0 



converges in L^(J7). 

Some unexpected results happen when we introduce a weight of the form 
f{Bf./n) in (|5.4p . In fact, a new critical value {H = ^) plays an important 
role. More precisely, consider the following sequence of random variables: 

^ n—l 

= 7^ E f{Bk/n)Hg{n^AB,/J. (5.5) 

Here, the integer g > 2 is fixed and the function / : M ^ M is supposed 
to satisfy some suitable regularity and growth conditions. In [13 HB], the 
following convergences as n — > oo are shown: 



• UH < ^, then 



„.«-iG„'il5'(l^P(.)(B.M,. (5.6) 



c.'^^' r nn.m.. (5.7) 

Jo 



where is a Brownian motion independent of B, and 



• If = 1 - ^, then 



Vlogn V 9! V 2gy V qj 

where is a Brownian motion independent of B. 
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• liH >1- ^, then 



Jo 

where Z^"^^ denotes the Hermite process of order q canonically con- 
structed from B (see |18j for the details). 

In addition, when q = 2 and H = ^, it was shown in [19] that G„ converges 
stably to a linear combination of the limits in (j5.7p and (j5.6p . (The proof of 
this last result follows an approach similar to the proof of our Theorem 13.11 
and allows to derive a change of variable formula for the fBm of Hurst index 
J, with a correction term that is an ordinary Ito integral with respect to a 
Brownian motion that is independent oi B.) But the convergence of G„ in 
the critical case H = q > 3, was open till now. 

In the present paper, we are going to show that Theorem 13.11 provides a 
proof of the following new result, valid for any integer q > 2 and any index 

G'„-n-i-''^i^5;/(<^)(S,/„)^^^^^aH,J f{Bs)dWs. (5.9) 
^' fc=0 "^^ 

(See Theorem 15.31 below for a precise statement.) Notice that (j5.9p provides 
a new proof of (j5.7p in the case H G (^i^) (without considering two 
different levels of discretization n ^ m, as in [18j). More importantly, in the 
critical case H = convergence (|5.9p yields: 

-^^Q- Jo Jo 

Hence, the understanding of the asymptotic behavior of the weighted Her- 
mite variations of the fBm is now complete (indeed, the case -ff = g > 3, 
was the only remaining case, as mentioned in the discussion above). 

The main idea of the proof of (j5.9p is a decomposition of the ran- 
dom variable Gn using equation (|2.6p . The term with r = is a mul- 
tiple Skorohod integral of order q and, by Theorem 15.21 below, it con- 
verges in law for any H £ (^^'|)- The term with r = q behaves as 

— n~2~'^^ ^"Z^-*! Ylk=o f'''^^ (^k/n)- The remaining terms {I < r < q — 1) 
converge to zero in L^{Q). 
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5.2 Some preliminaries on the fractional Brownian motion 

Before proving (|5.9p . we need some preliminaries on the Malliavin calculus 
associated with the fBm and some technical results (see \20\ Chapter 5]). 

In the following we assume H ^ (Oj |)- We denote by £ the set of step 
functions on [0,1]. Let S) be the Hilbert space defined as the closure of £ 
with respect to the scalar product 



+ t 



H 



t 



The mapping ljo,t] Bf can be extended to a linear isometry between the 
Hilbert space S) and the Gaussian space spanned by B. We denote this 
isometry hy (p ^ B{(j)). In this way {B((j)), (j) £ S)} is an isonormal Gaussian 

H-- 

space. (In fact, we know that the space coincides with Iq^ ^(L^[0, 1]), 
where 

H-h 1 

fix) 



is the left-sided Liouville fractional integral of order H — ^, 
Prom now on, we will make use of the notation 



see 



k/n 



'(fc+l)/ra ~ ^k/n — l(fc/n,(fc+l)/n]) 



for t £ [0, 1], n > 1, and A; = 0, . . . ,n - 1. Notice that Hq{n^ABk/J 
We need the following technical lemma. 



Lemma 5.1 Recall that H < ^. Let n > 1 and k = 



, n — 1 . We have 



(a) \E{Br{Bt - Bs))\ ^ (t - s)^^ for any r e [0,1] and ^ s < t ^ 1. 



(h) {et, dk/n)^ 



< n 



-2H 



for any t G [0, 1]. 



(d) For any integer q>2, 

n-l 



0(1) as n tends to infinity. 



E 

k=0 



0{n 



-2/^(g-l)^ 



as n tends to infinity. 

(5.10) 
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(e) Recall the definition i5.2\) of pn- We have 

{dj/n, dk/n)^-^ = n^'^^ pH{k- j). 
Consequently, for any integer q>l, we can write 

d^in)^ = 0{v}~'^'^^) as n tends to infinity. (5.11) 

k,j=Q 

Proof. We have 



E{Br{B,-Bs)) = l{r'^ + t'^-\t-r\'")-^{r'^ + s'"-\s-r\'^) 
= lit'''-s'^) + l{\s-r\^^-\t-rn. 

Using the inequahty {b"^^ — a"^^] ^ |6 — for any a, 6 G [0, 1], we deduce 
(a). Property (b) is an immediate consequence of (a). To show property (c) 
we use 

(^*' 9k/n)s, = ^ [(^ + l)'"" -k^^'-lk + l- nf ^ + \k- nf ^] . 
Property (d) follows from 

and 



dk/n) 



1 



1=1 



2H u2m 



< 



i=l ^ ^ 



Finally, property (e) follows from 



n-l 
fcj=0 



n-l 



fcj=0 reZ 
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5.3 An auxiliary convergence result 

From now on, we fix g > 2 and we make use of the following hypothesis on 
/ : M ^ M: 

(H) / belongs to C^'^ and, for any p > 2 and i = 0, . . . , 2q, 

E{snp |/»(i3i)n<oo. (5.12) 

te[o,i] 

Notice that a sufficient condition for (j5.12p to hold is that / satisfies 
an exponential growth condition of the form |/'-^'^^(a;)| ^ ke'^^^^^ for some 
constants c,k>0 and < p < 2. 

The aim of this section is to prove the following auxiliary convergence 
result. 

Theorem 5.2 Suppose H S ^^i^); '^''^d let f be a function satisfying 
Hypothesis (H). Consider the sequence of q-parameter step processes defined 
by 

n-1 

«n = n''^-i^/(i?,/„)af;„. (5.13) 
fc=0 

Then Un G DomJ'', and 5'^{un) converges stably to crH,q Jq f{Bs)dWs, where 
W is a Brownian motion independent of B, and (TH,q > is defined in \5. 

Proof. The fact that u„ belongs to Domd'^ is a consequence of the 
inclusion W^'^{f)®'i) C Dom^^ and hypothesis (H). We are now going to 
show that the sequence F„ = 5'^{un) satisfies the conditions of Theorem 13. 1[ 
We make use of the notation 

= {^k/n-,dj/n) ) A,i = (c?fe/„,5j/„)^ , (5-14) 

for /c, J = 0, . . . , n — 1 and n > 1. Also C will denote a generic constant. 

Step 1. Let us show first that -F„ is bounded in L^(0). Taking into 
account the continuity of the Skorohod integral from the space D'''^(i3®'^) 
into L'^{VL) (see (g^])), it suffices to show that u„ is bounded in D'?'2(^®'?). 
Actually we are going to show that n„ is bounded in J}^''p[^®^) for any 
integer k < 2q and any real number p > 2. Using the estimate ()5.1ip we 
obtain 

n-1 

\\un\\l^.=n^'"'-' V f{Bk/n)f{Bj/n)PL<C sup \f{Bt)\\ 
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Moreover for any integer k > 1, 



n-1 



j=0 



and we obtain in the same way 



n-1 



^0(9 + *) 



n 



1,3=0 

< C sup l/e^n^t)'- 



0<t<l 

Then the result fohows from hypothesis (H). 

Step 2. Let us show condition (i) of Theorem I3.1[ Fix some integers 
r, fci, . . . , kq-i > such that ki + 2k2 + . . . + {q-l)kq-i+r = q. Let h G Sj^"^. 
We claim that (n„, (DFn)'^^' ... O h)^^^ converges to 

zero in L^{Q). Suppose first that r > 1. Without loss of generality, we can 
assume that h has the form g (gi et, with g G S)'^^^^^\ Set = {DFn)'^^^ ® 
[D'i'^ Fn)®^"-^ ® g. Then we can write 

n-1 



As a consequence, 



fc=0 



n-1 



k=0 



Condition (c) of Lemma 15. II implies 



n-1 



A;=0 



< C. 



Hence, 



-nil fi(g>(cj-l) 



21 



On the other hand 

9-1 

llrFi l|2 _||„||2 TT ||n"^Z? Il^'^n 

\\^n\\^ig)(q~l) — ||y||^®(r-l) J_J_ W-L^ -f^nll fs®™ ; 

m=l 

and applying the generahzed Holder's inequality 

9-1 



9-1 

2 ^ ^ TT /'l? /'ll n»n Z7 ||2fcm(9-l)^ 

m=l 



9-1 

I I II T^mTp ||2fcm 

m=l 

By Meyer's inequalities ()2.9p . for any 1 < m < g — 1 and any p > 2, we 
obtain, using Step 1, that 

\\D FriWiPl^Q-sj®"^) ~ 11^ ^''(^n)llLP(r2;i5®'") 

— Il'^'^(''^n) IIb"*'? — ^ ll''^ri||]D)™+<j,p(f,«q) ^ C. 

Therefore, 

{\{Un,^n®et)^tSc\) < Cn^~^, 

which converges to zero as n tends to infinity because H < ^. 

Suppose now that r = 0. In this case, we have = {DFn)®^^ 
(L'''-ii?„)®fc,-i. Then 

»®;,$\ =(d.,n,DFnf' ■■■(d''j''-^\D'i-^Fn)^'''" . (5.15) 
From (j5.15p and (j5.13p we obtain 

n— 1 g— 1 ^ 

K,<fn)^«, = n-^^-^ n . (5.16) 

fc=0 m=l 

Notice that for any m=l,...,(7 — 1, the term ( c?*^,™", D'^Fn ) can be 
estimated by n~"^^ \\D"^Fn\\^(g,m- Then, taking into account that 



supS(p"^F„||P ) <oo 



for any p> 2, and that ^^^=1 "^^m = obtain for £^ ^n)sj«>q\) an 

estimate of the form Cy/n, which is unfortunately not satisfactory. For this 

reason, a finer analysis of the terms {d®J^,D^Fn) is required. 

\ J/" / f,®m 
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First we are going to apply formula (j2.1ip to compute the derivative 



E 

i=0 



m\ I q 



I I \ i 



n 



qH- 



i=0 



n-1 



_(55(r7i— j) 
''l/n 



I In 







(5.17) 



Set K'^ = (^dJJ^,D"'Fn)^^^, and recall the definition of akj and pkj 
from ()5.14p . From (j5.17p we obtain 

n-l 



1=0 ^ / V / 

m 

(5.18) 



i=0 



with 



1 I m 
i J \ i 



n-l 



5 (<?-«) 



«=0 



By Meyer inequalities (j2.9p we obtain, using also assumption (H), that, for 
any p>2, 



(/(— ')(i?;/„)9^'^. 



LP 



< c 



/(™-)(i?,/„)af(''-^) 



l/n 



ID)'J-».P(f)'8'J-i) 

(5.19) 



Using Lemma[0(b) and (e) we have a™r* < c'„-(™-*)2H ^nd 
(j^~2tH ^ Therefore, for any i > 1, we have 



< 



n-l 

II^^^'^'^ILp < Cn'"-^ J2 ^ C-^"^"^™^^'-^- (5.20) 

On the other hand, if i = 0, Lemma |5. II (c) and (j5.19p yield 
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Notice that the estimate for the L^(r2)-norm of <I>n™''-' in the case i = is 
worst than for i > 1. We will see later that, for p = 2, we can get a better 
estimate for $n 

Because J2m=i ^ra > 2, the number of factors in Hmii {^j/n-. D'^Fn)^^ 
is at least two. As a consequence, we can write 

n-l 
j=0 

for some fx, v (not necessarily distinct), where 

Qi = {^^^^f^^\n^)''"-^ n (*n'^y'"- (5-22) 

m=l 

Consider the decomposition 
where 

71—1 ( ii- y \ 

j=0 \i=0 k=0 / 

n-l 

j=0 

From (I5:22|) and the estimate H^'n'^ HiP < Cn'"^^ , for all p > 2 and 1 < 
m < q, we obtain 

\\Qi\\Lr < Cn^^^'^-^'^^'l (5.23) 
Then, from (lOOl) . (lOTT) and (lOHl) we obtain 

i=l k=l 

fj. V 
_^ ^ ^~l~2{ii+u)H+iH+2H _^ ^ ^-l-2(/i+i/)/f+fc//+2// j 

i=l fc=l 

which converges to zero as n tends to infinity, because /i, > 1 and H < ^. 
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For the term using again the estimates (j5.2ip and (j5.23p we get 



< Cn-2 



1+2H 



which converges to zero as n tends to infinity if < |. To handle the case 
i/ G [|, ^) we need more precise estimates for the L^(r2)-norm of $n'^'"'. We 
have, using formula (|2.12|) 



E 



2 / \ 2 

m 



n-l 



n 



2 / \ 2 n-l 
2gH-l Q] (m\ 2 



=0 



Ll'=0 



n 



2gH-i Q\ rn\ 2 



2 ra-1 



-io2i 

Plv 



«,/'=0 j=0 ^ ^ 



i=0 



If i > 1, then X]rz'=o A^l' — ^'^^ ™^ obtain an estimate of the form 
||-R^||^2 ^ CrC , where 

7 = ^ {2qH - 1 - 4i/i7 - A{q -i)H +1- AiH) = -qH - 2vH. 

For i = 0, then sup„ Yllv=o \ < c«, and we get 

7 = ^ (2gi7 - 1 - 2H{2v + 2q- 2)) = -qH - 2uH -^ + 2H. 

We have obtained the estimate 

K"'mL2 < Cn-'^H~2uH+2H-^ _ (5 24) 

Fix < a < 1. This choice is possible because < 1. We have, by 
Holder's inequality. 



n-l 



j=0 



25 



Using ([Oi]) . ([OT]) and (f5:23]l we obtain 

E {\Bn\) < Cn"' , (5.25) 

where 

7 = + ^ + [-2gF - 2(/i + + 4// - 1] a 

-//(g - /i - I/) + (1 - a)(-l - 2H{fi + + AH) 
= -^+AH-H{fi + u)-2aqH 

< -^ + 2H -2aqH <^-2aqH <0, 

because H < ^. Therefore E {\Bn\) converges to zero as n tends to infinity. 
Step 3. Let us show condition (ii). We have 

From (j5.18p we get 

9 / \ 2 n-l 
i=0 ^ ^ 1=0 

Therefore, we can make the decomposition 

{Un, D'^ En) f^igiq = An + Bn + C^, 

where 

n-l 

>1„ = n"i^-\\Y,fiy{Bi/n)f{B,/n), 
i,j=o 

g-1 / N 2 n-l 



i=l ^ ^ i,i=0 
n-l 



1,3=0 
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The term An converges to a nonnegative square integrable random variable. 
Indeed, 



n-l 



^« = ^Y.f^^k/n)f{B,/n){\k-j + ir + \k-j-l\^^-2\k-j\''y 
k,j=0 

. oo (n-l)A(n-l-p) 

= i^H E fiB,/MBu^,yn){\P+ir+\P-ir-2\pn' 

p=—oo j=OV-p 

which converges in L^{Q,) to 



\kez J 



Then, it suffices to show that the terms i?„ and C„ converge to zero in 

-2/o\ il„„j- ™ D ;i„ ,,„; f„„j- j-l,„j- V^™-1 



L (f2). For the term i?„ we can write, using the fact that Yl?j=o 



< 



g— 1 n— 1 

E E 



L2 



which converges to zero as n tends to infinity. Finally, for the term C„ we 
can write 



and i - 2qH + {2H - i) V < 0, because if - ^ < this is true due to 
\ - 2qH < 0, and if 2i7 - i > 0, then we get 2F(1 - g) < 0. This completes 
the proof of Theorem 15. 2[ ■ 

5.4 Proof of the stable convergence ( 15. 9j) 

As a consequence of Theorem we can derive the following result, which 
is nothing but (j5.9p : 
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Theorem 5.3 Suppose that f is a function satisfying Hypothesis (H). Let 
Gn be the sequence of random variables defined in \5. 5|) . Then, provided 
H G i), we have 

where W is a Brownian motion independent of B and crH,q > is defined 
by (E3j. 

Proof. We recall first that Hg [n^ (ABfc/„)) = ^n^" 5'> {d®JJ . Then, using 
yields 

r=0 ^ ^ 

where ak^k is defined in (j5.14p . As a consequence, 

^' r=OA;=0 ^ ^ 

1 '"^ 

^' r=l 

where m„ is defined in (j5.13|) . 

^" ^ ^ A;=0 

and 

^ n— 1 

^' k=0 

The proof will be done in two steps. 

Step 1 We first show that if G (O, i) , and r = 1, . . . , g - 1, 5'^-'' {v'^n^) 
converges to zero in L^(J7) as n tends to infinity. It suffices to show that Vn 
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converges to zero in the norm of the space D'' ^''^(^S^'^^i ^)). For < m < 
q — r, we can write, using the notation ^ defined by (j5.14p . 



E 



^®((j-r + m) 



1 

q\ \r 

n-1 

X 

k,l=0 



n 



2qH-l 



J2 E(f(^-+"'\B,/^)f^^+^\B, 



l/n, 



I 

< (j^2qH-l^-2H{2r-2+m+q-r) 



which converges to zero as n tends to infinity. 

Step 2 To complete the proof it suffices to check that 



n-1 



fc=0 



converges to zero in L?{^1) as n tends to infinity. This follows from (15.10p 
and the estimates 



n-1 



n-1 



fe=0 

n-1 



2<iq\ 



k=0 



L2 



k=0 



1 



a 



k,k 



Notice that -qH + 2H - \ <Q. The proof is now complete. 
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